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The high-density electron gas in a strong magnetic field B and at zero temperature is investigated. 
The quantum strong-field limit is considered in which only the lowest Landau level is occupied. It 
is shown that the perturbation series of the ground-state energy can be represented in analogy to 
the Gell-Mann Brueckner expression of the ground-state energy of the field-free electron gas. The 
role of the expansion parameter is taken by tb = (2/3-7T 2 ) (B/m 2 ) (hrs/e) 3 instead of the field-free 
Gell-Mann Brueckner parameter r s . The perturbation series is given exactly up to o(r_g) for the 
case of a small filling factor for the lowest Landau level. 
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O : I. INTRODUCTION 

o 

The calculation of the correlation energy of a field-free and fully degenerate electron gas has a long history. Since 
the pioneering work of Gell-Mann and Brueckner [jlj and Wigner who derived analytic results for the high-density 
and the low-density limit, respectively, various interpolation formulas between these two limits have been established. 
For a free electron gas the ground-state energy is a function of the interelectron spacing r s only, which is related to 
the particle density n by n^ 1 — (Aitr^a\l?>), where as is the Bohr radius. Especially for a high-density electron gas 
the ground-state energy takes the form Tit] 
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O' e g = — '— + 0.0622 ln(r s ) -0.094. (1) 
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The leading term is the kinetic energy. The next term is the first-order exchange energy, while the remaining terms 
. in this series are called the correlation energy. 

The purpose of our work is to find an analogous expansion of the ground-state energy for a system of a large number 
of electrons moving in a fixed uniform distribution of positive charge and in an external uniform magnetic field. In 
I ■ particular, we focus on the strong-field limit, where all the electrons are in the lowest Landau eigenstate and the spins 
are all aligned antiparallel to the magnetic field. The system is assumed to be at zero temperature. In our analysis 
we essentially follow the calculation originally developed by Gell-Mann and Brueckner. A homogeneous magnetic 
field B modifies the energy spectrum of a charged particle and the Fermi energy (i.e., the chemical potential of the 
noninteracting magnetized electron gas) is no longer a function of r s alone, but also of the magnetic field. Therefore, 
it is convenient to introduce the new parameters 
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where lu c = eB/m is the cyclotron frequency, a = as/Is is the ratio of the Bohr radius, and the magnetic length 
Ib = y«7 {gB). At strong magnetic fields tb takes the role of the expansion parameter. The second parameter t may 
be regarded as a filling parameter. It must satisfy the condition t < 1 in order to describe a quantum system in the 
strong- field limit. In both equations we have made use of the the relation between the Fermi wave vector kp and the 
particle density in the strong- field limit, where the system behaves essentially as a one-dimensional electron gas 
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— kp 



dk z , k F = 2n 2 l B n. (3) 



Now the kinetic energy of the noninteracting degenerate electron gas per particle in rydbergs becomes 

I , ... , . 0.0337 

eun = ,-_,_, 2 / ^— dk z —— = —2-2- « —2—. (4) 
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Thus the kinetic energy is a function of the dimensionless parameter tb only. In the limit rs < 1, the kinetic 
energy will give the dominant contribution to the ground-state energy, while the Coloumb interaction acts as a small 
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perturbation to the motion of the electrons. We investigate here this physically interesting high-density regime. We 
show that in this case the ground-state energy takes the form 



e g = — 2 — H h B{t) ln(r_e) + C{t) + terms that vanish as re — > . (5) 

The second term in this series comes from the first-order exchange energy. Its calculation was given by Danz and 
Glasser and will be reviewed in Sec. II. The terms in the energy beyond the Hartree-Fock approximation can be 
obtained by a formal summation of all ring diagrams. In addition to that, we must include the second-order exchange 
diagram that contributes to C(t). Depending on the values of the magnetic field and the density the filling parameter 
t may vary within the range < t < 1. An analytical calculation for the coefficients A(t), B{t) and C{t) is too 
difficult as it involves a troublesome summation over Landau indices. Therefore, we give explicit expressions for the 
contributions to the constants A(t), B(t), and C(t) and find analytic results for their asymptotic behavior as t goes to 
zero. The validity domain of this expansion in the r s -t plane is schematically shown in Fig. 1, in which we have made 
use of the relationship tb ~ r s /ts to indicate the high-density region (tb <C 1). The regime of a coupled electron 
gas (r_e > 1) was investigated by Fushiki et al H using the Thomas-Fermi type statistical model. At sufficiently low 
densities the electron gas should form a Wigner crystal. Following Kleppmann and Elliott [Q (see also Ref. ||), a 
strong magnetic field would increase the density at which crystallization occurs. The results obtained in this work 
may be used to establish interpolation formulas between the various limiting results. 
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FIG. 1. Scetch of the validity domain of the various expansions in the r s —t plane: (a) high-density region [the cross-hatched 
area at the bottom shows the validity domain of the analytic results for the asymptotic behavior t — > of the coefficients A{t), 
B{t) and C {€).], (b) intermediate-density region, and (c) low-density region. 



The strong field limit has also been investigated by Horing et al. [|6| and Isihara and Tsai |J7|. We have partially 
adopted the notation of Horing et al. and additionally we have introduced the parameter t to clarify the structure 
of the perturbation expansion of the ground-state energy Eq. (^) . This was not explicitly given in Ref. Q . Horing et 
al. || carefully analyzed the spectrum of the plasma oscillations within the random-phase approximation to calculate 
the correlation energy. It is found to be proportional to ln(re). We show that the same result can be obtained by 
expanding the polarization function in powers of the momentum transfer, which gives a better understanding of the 
structure of the ground-state energy expansion and simplifies the calculation. Furthermore, this enables us to find 
the constant C(t) of this series. This will be discussed in Sec. III. 

The calculations of Isihara and Tsai H are performed in the grand canonical ensemble. Within our notations their 
ground-state energy agrees with that of Horing et al. ||] up to the order A(t)/rB, but differs in the constant B(t) by 
a factor 2. We confirm in our calculation the result of Horing et al. |(|. Isihara and Tsai @ also give contributions to 
the term C(t); however, within this order their calculation is not complete. 

Our work is motivated by the observation that all of these conditions [degeneracy, weak coupling {tb *C 1) and 
small filling factor (t <C 1)] may be realized on the surface of some strongly magnetized neutron stars ||. The 
outermost layers of these objects are likely to form a degenerate hydrogen plasma at magnetic field strengths of I0 9 T, 
at temperatures of lO 5 ^ and at densities of 3 x 10 34 m~ 3 . Then the expansion parameter and the filling factor are 
found to be tb ~ 0.01 and t f» 0.05. On the other hand, at higher temperatures (T > 10 5 -ftT) and lower densities 
(n < 10 30 m~ 3 ) the magnetized plasma (B = 10 9 T) is nondegenerate and weakly coupled. This situation has been 
recently considered in Refs. and p0| . 

The high-density strong-field limit may also be achieved in laboratory plasmas, such as semiconductors or laser- 
induced plasmas. In particular, one may consider semiconductors with small effective mass m* — 0.01m e and carrier 
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densities of 5 x 10 23 m -3 in strong magnetic fields of 40T, such as indium antimonide. At these values of B and n 
the parameters tb ~ 0.02 and t « 0.04 lie within our approximation scheme. Throughout this paper the energies are 
given per particle and in rydbergs. 



II. EXCHANGE ENERGY 



In this section we shall briefly review the calculation of the first-order exchange energy, which is diagrammatically 
represented in Fig. 2. 



FIG. 2. First order exchange diagram. 

This contribution has been examined by Danz and Glasser . Following the rules of calculating Feynman diagrams 
we have 



(-hf = 
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J d ld2V x (12)G<(l2)G<(21). 
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The trace is over the spin variable o z = (—1, +1). Throughout this work we shall use the closed form of the Green's 
function G<(12) in space-time representation 1 = (ri,ti) as given by Horing jl2), which has also been exploited in 
Refs. [§ and |ftl[ 



du J - /o(w)] 
ifo(u>) 



-iuiT 



dT' 



(2tt 



cos(^I* 



exp 



.pI+pI . 

-i tan 

mtOc 



(7) 



Here we have introduced the Fermi-Dirac distribution function fo(u>), the Bohr magneton its = eh/ (2m) and a unitary 
phase factor C(r, r'). After inserting expression (Q) into Eq.(|]) and considering the quantum strong-field limit only, 
i.e., huj c > ep, one can perform all elementary integrals to obtain the Hartree-Fock exchange energy 
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In the limit t — > the exchange energy takes the simple form 
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e HF (t -> 0) = ^(A-C-2 ln(2) - ln(t)) 
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where C is Euler's constant C w 0.5772. This result was already obtained by Danz and Glasser JTTJ . One can easily 
identify the coefficient Ait) from Eq.(||) and its asymptotic limit t — > from Eq.Q. 



III. CALCULATION OF THE CORRELATION ENERGY 



In second-order perturbation theory, one expects a higher-order contribution in the interaction parameter e 2 or 
t\b, that is, a constant independent of re- Like in the zero-magnetic- field counterpart, this contribution diverges 
logarithmically at small momentum transfers. As was shown by Gell-Mann and Brueckner, this difficulty can be 
overcome by summing ring diagrams up to infinite order (Fig. 3), which is known as the random-phase approximation. 
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FIG. 3. Ring approximation for the correlation energy. 

Given the polarization function II(p, u)) within the random- phase approximation, the correlation energy may be 
written as 
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Here and in the remainder of this paper we have expressed all momenta in terms of the Fermi wave vector k,F- 
For a further detailed analysis, the polarisation function II(p, u) must be evaluated. This quantity is determined 
solely by a collisionless, i.e., noninteracting, Fermi system in a magnetic field and may be calculated in a space-time 
representation by 



11(12) = ^?r 1 Tlv {G a {\2)G a {2\)) . 
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By using Eq. ([?]), the polarization function was studied in detail by Horing Q. He obtained useful analytical results 
for the nondegenerate as well as for the degenerate region. In particular he derived an exact expression for the 
polarization function in the quantum strong-field limit [Ref. |13], p. 61, Eq.(A.III.6)), which is given by 
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Here we have introduced our notation [Eq.(^)]. We have also made use of Eq.(||), which relates the Fermi wave vector 
hp to the expansion parameter tb- Equations ( pj] ) and ( [l3|) represent the basis of our further calculations. From 
the ring approximation e r we may find the contribution to B(t) and C(t) by expanding the polarization function 
in powers of the momentum transfer p. Only the lowest-order term independent of p must be considered, while 
higher-order terms proportional to p 2 will give higher powers of tb in the final expression of the ground-state energy. 
An additional contribution to C(t) arises in the exchange interaction of the order e 4 . This energy may be represented 
diagrammatically as in Fig. 4. 




FIG. 4. Second-order exchange diagram. 

Higher-order exchange terms lead to contributions in higher powers of tb and need not be considered in our 
calculation. 



A. ring approximation 



We may start our calculation with the investigation of the ring approximation. Following the original work of 
Gell-Mann and Brueckner, we first establish an expansion of the polarization function in powers of the transfer 
momentum. For this we introduce the dimensionless frequency u — > (rruo) / \Tikp\p z \) . Now we fix this quantity and 
let the momentum transfer approach zero to obtain the limiting form 
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In general the functions R\, R2, etc., will depend on the direction of the transferred momentum with respect to the 
magnetic field that is taken into account by the additional argument cos(#) in R2. However, as may be easily shown, 
Ri is independent of this angle and takes the simple form 



Ri{u) 



1 + u 
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Now we perform the integral over p in Eq. ( |ll[ ) from to 1 and correct to order r B exactly. This may be achieved by 
introducing S(t), which restores the second-order pertubation term to its exact value. Thus we have 
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This is a finite expression since the logarithmic divergence cancels. Substituting Eq. (|l3| ) into this equation, one 
obtains the general form of this S(t) as given in Appendix A [Eq.(Al)]. There we have derived an asymptotic 
expansion that is found to be 
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where we have introduced the Rieman zeta function £ with £(3) ~ 1.202057. The integration over the momentum in 
Eq. ( |l6| ) is readily carried out with the result 
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where we have dropped all terms that vanish as r B — > 0. The remaining integral may be performed exactly and we 
find for the logarithmic contribution, i.e., for B(t), the result 
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Note that B{t) is exact for arbitrary filling factor. This result is in agreement with that of Horing et al || . Furthermore, 
we are able to determine the next term in the perturbation series. Again, we can perform the integral (|l9| ) and take 
into account the contributions coming from S(t) and e e x(t) to obtain 



C(t) = 
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We may proceed with the calculation of the second-order exchange term in the next subsection. 



B. second-order exchange term 



This contribution may be written in space-time representation as 
1 „ f 1 d\ 



Eex{t) = 2^d Tr,7 l T I dl I d2 J d3 / d41/A ( 14 ) vrA ( 23 ) G -( 12 ) G -( 23 ) G -( 34 ) G -( 41 )- ( 22 ) 

Now we may use Eq.(Q) and employ the expansion of the Green's function in terms of the Laguerre polynomials L s (x) 
[Ref. Q, p.96, Eq.(2.3)]. Then the V and oj integrations are trivial and in the limit Huj c > cf this contribution can 
be written as 
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-6(s - l)0(s' - 1)L S _! [4t(p, - q p ) 2 ] L«,_i [4t(q p + k p ) 2 ] ) . (23) 



Here p p denotes the two-dimensional vector p p = (p x ,p y ). This expression is the basis for the caculation of C(t) at 
arbitrary filling factor t < 1. For small values of t we may seek an asymptotic exact result for e ex (t) . We find that in 
this limit only the state s — s' — contributes to the exchange energy, while the other terms (s > 1 or s' > 1) vanish 
as t goes to zero. Thus we have, for 
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The detailed evaluation of this integral in the asymptotic limit is given in Appendix B. It shows an interesting 
logarithmic dependence on the filling factor, which is given by 
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IV. CONCLUSIONS 



In this paper the ground-state energy of a degenerate electron gas in a strong magnetic field has been investigated. 
The electron gas in the strong-field limit is characterized by two parameters: the dimensionless inverse Fermi wave 
vector tb and the filling factor of the lowest Landau level t [see Eq.(||)]. We have shown that in the case of a dense 
electron gas <^C 1 the internal energy per particle (measured in rydbergs) may be expressed in the form of Eq. dq). 
The constant A(t) is fully determined by the Hartree-Fock exchange energy and may be identified by using Eq.(g|). 
The next term in this expansion comes from the random-phase approximation correlation energy. The constant B(t) 
of the In(r^) contribution is given by Eq. (|2C|) . For the calculation of the constant C(t) one needs to include the 



second-order exchange energy. Given 5(t) and e ex (t) by Eq.(Al) and Eq.(|23|), respectively, one may calculate the 
constant C{t) using Eq.(|T]) at any t. 

For the case of very strong magnetic fields with a small filling factor t (but sufficiently small tb) we have given 
analytic expressions for the constants A{t) by Eq.(^|) and B{t) by Eq.(pO|). Our results for A(t — > 0) and B(t — > 0) 
coincide with earlier calculations of Danz and Glasser fyl and of Horing et al. H . The calculation of the constant 



C(t — » 0) is, to our knowledge, given for the first time. We can use the asymptotic result for S(t) [Eq, 18 and for 
e ex (t) [Eq^q] to get an expression for small values of t: 

c(i ^ ) = -^-?if±i + ^fl + (-5 + | + M2) + ?£2) 

lon^t 16ir z t 6 \ 6 3 7r z 

One should emphasize that Eq.(||) with the asymptotic results for the constants [Eqs.(|9|),(p0|), and (p6|)] is an asymp- 
totic expansion of the ground-state energy of an electron gas in the limit re — * and then t — > with a possibly poor 
convergence. Including the results of our earlier work where the thermodynamic perturbation expansion in the 
Boltzmann limit has been performed, we have now found exact results of the thermodynamic functions in the limiting 
cases of a degenerate and nondegenerate weakly coupled electron gas. A third known limiting case is the region of 
Wigner cristallization (see Refs. Q and |J and references therein). One may now try to connect the known limiting 
results by a Pade approximation in order to describe the intermediate region as indicated in Fig. 1. A survey of 
several Pade approximations for the case of the zero-field electron gas was recently given by Stolzmann and Blocker 
p3| . A Pade approximation for the thermodynamic functions of the electron gas in the quantum strong- field limit 
will be given elsewhere. 
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APPENDIX A: CALCULATION OF S 



This contribution may be written in terms of the polarizationfunction ( |13| ) as 
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We are interested in small values of t, therefore we can restrict the calculation to n = m = as n or m > 1 give rise 
to higher-order terms in t. Note that the terms n or m > 1 are convergent expressions even without subtracting the 
logarith mic divergent part, while the lowest-order term n = m = would show singular behavior. The ui integration 
in Eq. (Al) may be performed exactly. The result may be split into small- and large- wave-number contributions 
5{t) = S 1 {t) + 8 2 (t) with 
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First we focus on the calculation of S\(t). The z integration may be expressed in terms of the exponential integral 
Ei(-x) 
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Now we can expand the logarithm and subtract the divergent part from the lowest-order contributions. Thus we have 
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Again, we can expand the exponential and the exponential integral and retain only contributions that do not vanish 
as t — > 0, 
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This may be simplified by expressing the sums in terms of the ip function and observing that they give elementary 
constants for integer arguments. Therefore, the final result may be written as 
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We have also introduced the Rieman zeta function C,(x). Following the same steps as discussed before, we obtain for 
the large momentum part 
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The transformation p 2 — > the expansion of the exponential and an integration by parts yield 
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By expanding the exponential integral the relevant contributions for the asymptotics are identified to come from the 
n — term. Then we find 
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Performing the sums exactly the result reads as 
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One can combine the results for d±(t — > 0) [Eq.( |A7D ] and (^(f — > 0) [Eq.( All] )] to find the asymptotic result for 
<5(t -> 0) [Eq.©]. 



APPENDIX B: CALCULATION OF THE SECOND ORDER EXCHANGE TERM 



After the substitution k' = (p — k — q) Eq. (K4) may be rewritten as 
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The Gaussian integrals can be done immediately. The analysis of the conditions on the region of integration gives 
then 
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Performing the p z and q z integration we find 
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We split up the integral into two contributions e ex (t) = Ii (t) + I 2 (t) with 



(B3) 



If f e~~ 2q e e~ 2kp 



dk z In 1 



k 2 t 



l„|l + ^]l„| i 



In 1 



^ ^ - In ( 1 + 

(2-k z ) 2 t) { + k 2 t 



(2 - k z )H 



(B4) 



and 

'■<" - -a? W^ 2 ? jf * h + i) { 21 ° + Mi) - 21 ° (' + ra)} ■ (B5) 

First we concentrate on Ii(t), where after the substitutions k p /k 2 — > etc., the k z integration is again of the 
Gaussian form. After the transformations k = q p k p /{q p + k p ) and q — k p J (q p + k p ) we have 

hit) = L ( dq r dk 3 C ^ = 1 = In f 1 + In fl + - ^ ^ . (B6) 

Now we investigate this integral in the limit of small t and find for the asymptotic behavior a logarithmic divergence. 
Noting that 
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—r ln(t) / / rifc In - = (3 ln(2) - 2 + C) ln(t) , 

/ /°° rlh e 4fc - 1 _ h ( M i n C fc ^ ) = f 4 - 2C + — + — + 3Cln(2) 

-61n(2) + ^ln 2 (2)-^) (B7) 

and neglecting all terms that vanish as * — > 0, this contribution becomes 

/i r;2 2 q \ 

7i(* -» 0) = - ln 2 (i) + (3 ln(2) - 2 + C) ln(t) + 4-2C+ — + — + 3C ln(2) - 6 ln(2) + - ln 2 (2) . (B8) 
\2 2 6 2 J 

We now find an appropriate approximation scheme for the calculation of I2 ■ We may perform the identical k and q 
integration [see Ref. JTj], p. 530, Eq.(6)] to obtain the result 

Hi) = -±\ f dk z L\2ktt) + 2C ln(2. 2 *) + C 2 + £ - 2 £ ^l±m^& ( 2 ^)« ] 

" L n=l ) 

x {2 (fc 2 «- kl) ~ 2 (/c 2 «- (2 - k z ) 2 )} . (B9) 
We find this may be exactly rewritten as 

I 2 {t) = CI ln 2 (t) + C2 ln(<) + C3 + (C4(n) ln 2 (t) + C5(n) ln(*) + C6(n)) *™ (BIO) 

n=l 

with the constant coefficients 

r-2 



1 

~ ~3' 

C2 = JL jf dfc, ^2 ln(2A: 2 ) (2C + ln(2fc 2 )) (ln(2[2 - k z ] 2 ) - ln(2fc 2 )) + 2 ln(2fc 2 ) (Bll) 

x (ln 2 (2[2 - k z ] 2 ) + 2C ln(2[2 - ^] 2 ) - ln 2 (2fc 2 ) - 2C ln(2/c 2 )) j 

= |-|C + |C(3)-2 M 2) > 
1 f 2 

C3 = / dfc 2 ln(2fc 2 ) (2C + ln(2fc 2 )) (ln 2 (2[2 - fc z ] 2 ) + 2C ln(2[2 - A:^ 2 ) - ln 2 (2fc 2 ) - 2C ln(2/c 2 )) 

, ^ oom.^ C2 8 2 _ 2 4^ , C(3) 



-31n 2 (2) + (4 - 2C) ln(2) - — - - - -tt 2 + -C + ^ (241n(2) + 8C - 16) . 

3 3 90 3 7r z 



(B12) 
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We neglect higher-order terms in this calculation. Thus we get the asymptotic expansion 

hit 0) = -\ ln 2 (t) + - He + ^C(3) - 2 ln(2)) ln(t) - 3 ln 2 (2) + (4 - 2C) ln(2) - 

-\ -^ 2 + t C + ^ (241n(2) + 8C ~ 16) ■ (B13) 

Summing I\ and I2 , we have the final expression for the second-order exchange energy at small values of the filling 
parameter t given in Eq. (p5|) . 
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